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Aim 

To derive the macroscopic properties of paramagnetic salts and 
the 3D simple harmonic oscillator. 

Objectives 

1. To describe the basic physics of a paramagnetic salt. 

2. To use the partition function to calculate thermodynamic 
quantities. 

3. To derive the magnetisation of a paramagnetic salt. 

4. To explain the use of the 3D simple harmonic oscillator to 
model thermal vibration. 

5. To explain the equipartition theorem. 
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3.1 Paramagnetic salts 
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Paramagnets 



We have derived the Boltzmann distribution in theory. It has 
also be widely observed in experiments. 

One good example is the measurement on paramagnetism 
carried out by Warren E. Henry in 1952, at the Naval Research 
Laboratory in the United States. 




He placed salts of chromium, iron or gandolium in liquid Helium 
and measured their magnetisation in a magnetic field. 

Then he calculated the magnetisation with a formula that is 
derived using Boltzmann distribution, and found that the 
results agree very well. 
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Warren's results, published in Physical Review 



On the left is the experimental setup used by Warren. It 
contains liquid helium to cool down the paramagnetic sample 
near the bottom. 




highest value of HIT. 
Fig. 1. Schematic diagram of assembly of meLal Dewar for liquid 
helium and sample; displacement lift is shown at top. 



On the right are the measurements and calculations which 
agree well. 
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A Spin-1/2 Paramagnet 



We shall learn how to do this type of calculation. 

As with Warren's experiments, we consider a solid, 
paramagnetic salt. Such a salt contains ions that have 
magnetic moments. The magnetic moment in each ion comes 
from the electrons in the ion. 

Consider a salt with an angular momentum of J for each 
magnetic ion. 
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A Spin-1/2 Paramagnet 



We shall consider the simple case where J = S = 1/2. Then 
the magnetic dipole moment the ion is given by 

where \±b is the Bohr magneton . 

The Bohr magneton is a unit of measure for the magnetic 
moment of the electron. It has the value: 

^ B = g.274 x 10~ 24 J T" 1 . 

Suppose that a magnetic field B is applied to the salt. Since 
the ion has a magnetic moment, it would have a magnetic 
potential energy given by 

e = — /i.B = —fij^B or + /i/^ii 
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A Spin-1/2 Paramagnet 



The sign of the magnetic potential energy 

e = — /i.B = —fij^B or + njjB 

would depend on whether the ion is spin up or spin down, as 
shown in the figure to the left. 



e x = -flB 



s 2 = juB 

o 



Energy 



As = 2juB 



The magnetic field thus splits the energy into 2 levels. We shall 
call this a spin-1/2 paramagnet. 
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A Spin-1/2 Paramagnet 



In the spin-1/2 paramagnet, some ions are spin up, and others 
are spin down. 

a © © © (D 
= © © © © 
- © © © © 
© © © © 

Suppose there are TV magnetic ions, with a fraction p spin up 
and a fraction q spin down. 

The magnetisation is then defined by the resultant magnetic 
moment: 

M = pNfji B - qNfji B = Nfji B (p - q) 

Our result on the most probable macrostate would suggest 
that they follow the Boltzmann distribution 

rii = Aexp(-Si/k B T). 
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A Spin-1/2 Paramagnet 



Since 

e = — /i.B = —fi B B or + /i/jii 
the Boltzmann distribution is 

n\ = Aexp(n B B /k B T) and n 2 = Aexp(—n B B/k B T) 
To simplify the expression, define 

x = n B B /k B T. 
Then rewrite the Boltzmann distribution as 

ni = Ae x and n 2 = Ae~ x 

If the magnetic dipoles follow this distribution, the fractions 
would be given by 

ni n 2 
p — — and q — — 

N N 
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A Spin-1/2 Paramagnet 



The fractions are 



ni n 2 
p = — and q = 



N N 
Using the Boltzmann distribution, the total number is 

N — n\ + 712 

= A(e x + e~ x ). 

The fraction then becomes 

e x 

V = 



e x + e~ x 

and 

e -x 
q = 

e X _|_ e X 
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A Spin-1/2 Paramagnet 



The magnetisation is defined by 

M = Nn B (p-q) 
Substituting our results on p and p gives 



rf* *y* 

e x e x 



M = Na B , 

I J—S I rf* I rf* rf* I rf* 

e x _|_ e x e x _|_ e X 

= TV/^tanhx 



Substituting 

x = i_i B B/k B T, 

we get 



M = 7Vu B tanh f ) 



This is the equation for magnetisation that we wanted to find 



PHYS393/Hock 



12 



2009-2010 



Paramagnets with Angular Momentum J 



For paramagnets with momentum J that is more than 1/2, 
there would be more than 2 energy levels. The same principles 
can be used to find the magnetisation, and the result is 



M = N gjHB 



(J + 1/2) coth[( J + 1/2)C] - S£ th ( C / 2 ) 



2 

where C = gji^B^ /k B T and gj is the Lande g factor. 
This equation gives good agreement with experiments. 




Vr* 10"' SWISS /DEC 
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Partition Function 



We now look at how to apply the Boltzmann distribution to 
find the heat capacity and entropy of a paramagnetic salt. 

The heat capacity behaves in a distinctive way at very low 
temperatures. The entropy is important for low temperatures 
cooling techniques. We shall see more of these in the second 
part of these lectures. 

In order to calculate them, it would be useful to develop the 
method of the partition function. 
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Partition Function 



We shall first develop the method of partition function. This is 
essentially a series of formulae that will make it easy to 
calculate various quantities in thermodynamics from the 
Boltzmann distribution, 

rii = Aexp(-Si/k B T) 
To see how this works, we start by defining the 

Z = Y / exjp(-e i /k B T) 

i 

This is just a sum of the Boltzmann factor over all energy 
levels of the system that we are interested in. In this case, that 
would be the paramagnetic salt. 
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Partition Function 



Next, look at the total energy 

i i 

Compare this with the definition of the partition function Z 
above. The terms that are summed look similar except for the 
factor of Si. 

The idea now is to make a formula that can allow us to 
calculate the energy U from the partition function Z. Try 
differentiating the partition function 

Z = ^exp(-e < /fc B T) 

i 

with respect to temperature T. This gives 



= ^2L £ i ex P(- £ ^) 



dT 
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Partition Function 



Rearrange the derivative 



= ^2L £ i ex P(- E i/W 



to make the summation the subject: 



^ £i exp(- £i /fc B T) = fe B T 

We now have the summation part of the energy 

U = AJ2^ exp(-s,//e B T) 

i 

in terms of the partition function Z. 

Next, we want to express the factor A in terms of the partition 
function as well. 
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Partition Function 



To relate A to the partition function, consider the total number 
of particles: 

N = = Aj2exp(-Si/k B T) 

i i 

Notice that the summation is just the partition function Z 
defined above. Therefore 

N = AZ or A = N/Z. 

We can substitute this and the previous result 

dZ 

$>;exp(- £i /fc B T) = k B T 2 — 

i 

into the formula for the energy 

U = Aj2si exp(si/k B T). 

i 

This is the result: 

u = Nk B T 2 dZ 
Z dT 
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Partition Function 



Notice that in this result 

u = Nk B T 2 dZ 
Z dT 

the total energy is expressed in terms of the partition function. 

The reason for doing this is to make it easy to calculate U. In 
the paramagnetic salts and other systems we shall look at, it is 
often possible to obtain a fairly simple formula for the partition 
function, which can then be substituted into the above formula 
to find the energy U. 

We shall see that it is possible to find similar formulae for other 
quantities, such as entropy, magnetisation and Helmoholtz free 
energy, in terms of the partition function. 

The above formula for the energy can also be written as 

9 d\nZ 

U = NkuT 

a dT 

You can check this by differentiating the InZ. 
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Partition Function 



Next, consider the entropy 

S = k B \nn 

We want to turn this into a formula in terms of the partition 
function as well. 

In the macrostate (ni,n2,...), the number of microstates is 
given by 

Q = . 

Ui rii\ 

Applying Stirling's theorem gives: 

In Q = TV In TV— TV— (ni In ni— n{) — (n<2 In n^—n^) — ... — {n i In n^— rij) — ... 

Substituting into the entropy formula above and collecting the 
terms into a sum gives: 

S = k B (N\n N - J^njn 

i 
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Partition Function 



To express this 

S = k B (N\n TV - In n i) 

i 

in terms of the partition function, start by substituting the 
Boltzmann distribution: 

rii = Aexp(-Si/k B T) 

This gives 

S = k B [N\nN -^2ni0nA-£i/k B T)] 

i 

1 

= Nk B \nN - k B \nAY,n i + -J2 n i £ i 

i i 

The sums in the 2nd and 3rd terms are respectively the total 
number of particles TV and energy U. Therefore 

S = Nk B In TV — Nk B In A + ^. 
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Partition Function 



We have obtained 

S = Nk B In TV — Nk B In A + ^. 

Recall that we already have the formulae for TV and U in terms 
of the partition function: 

N = AZ and U = Nk B T 2 ^^ 

Substituting these, we get 

S = Nk B \nAZ -k B N\nA + Nk B T dlnZ 



dT 

= Nk B \nZ + Nk B T^^ 
This is the formula for S in terms of the partition function Z. 
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Partition Function 



We would need one more formula, for the Helmholtz free 
energy: 

F — U — TS. 

We shall first review the results that we are going to use. From 
thermodynamics, we know that 

TdS = dU + pdV. 

where TdS is the heat absorbed by a system, dU the increase in 
internal energy, and pdV the work done by expanding against 
an external pressure. 
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Partition Function 



Taking the differential of the Helmholtz free energy: 

dF = dU - d(TS) 
and substituting for dU from the previous equation gives: 

dF = TdS - pdV - d(TS). 

Simplifying, we get 

dF = -SdT - pdV. 

This equation makes it convenient to think of F as a function 
of T and V. We know from partial differentiation that the 
differential of the function F(T,V) is 



Comparing this with the previous equation, we see for example 




that 
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Partition Function 



We have found that the relation 




= -P- 



So if we can express the Helmholtz free energy F in terms of 
the partition function, we can use it to find the pressure p. This 
example would be useful later on in the lectures on ideal gas. 

We have previously obtained the formulae for U and S in terms 
of the partition function: 



Substituting these into the definition for Helmholtz free energy: 



U = Nk B T 



2 <9lnZ 



dT 



S = Nk B \nZ + Nk B T 



d\nZ 



dT 



F = U -TS, 



we obtain 



F = -Nk B T\nZ. 
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Formulae related to the partition function 



Boltzmann distribution 

rii = A exp ( — 

Partition function 

Z = £exp(- 

i V 

Total number of particles 

N = AZ 

Total energy 

9 <9lnZ 

U = NkuT 

B 8T 

Entropy 

d\nZ 



S = Nk B \nZ + Nk B T 

ol 

Helmholtz free energy 

F = —Nk B T In Z 
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Thermodynamics of the spin-1/2 paramagnet 



We can now find the internal energy spin-1/2 paramagnet. 

For the spin-1/2 salt, the energy levels are ~hb b ancl +VbB 
So the partition function is 

-Eexp(-^)=e XP (g|) + e X p(-g) 
This can be written as 

Z = 2cosh (^]. 

\k B TJ 

To find the internal energy, the formula is 

U = Nk B T 2 d\nZ/dT 
Substituting the above expression for the partition function 



gives 
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Behaviour of the internal energy 



The behaviour of the internal energy at high or low 
temperatures gives insight into the physics of the paramagnet. 



U = -Nu B Btanh [ f^.) 

\k B T) 




When U is plotted against T, the graph goes to -Nii B B for 
small T, and zero for large T. 



Can we explain this physically? 



PHYS393/Hock 



28 



2009-2010 



At low temperature, most atoms would be in the ground state 
- the lower energy level —\± B B. 

+HbB 

The total energy would be TV times this, which is —Nn B B. 

At high temperature, the difference between the energy levels 
2f_i B B would be small compared with the k B T in exp(—e/k B T). 

+H B B — | 

-mB — [ 

So the levels would be nearly equally populated: half at —\± B B, 
half at +fi B B. The total would go to zero. 
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The heat capacity 



We have previously obtained the energy 



U = -Nu B Btanh [ f^.) 

\k B T) 



We can obtain the heat capacity by differentiating with respect 
to temperature T: 

c = du 



dT 

This gives: 
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The heat capacity 



For high temperature - i.e. when k B T is large compared with 
\i B B - the argument in the hyperbolic secant function becomes 
small: 

k B T 

and so the hyperbolic secant tends to 1: 



sech 



\k B Tj 



So at high temperature the heat capacity becomes 

I II. D F> \ ~ 

C Nk B 



\k B TJ 



which is proportional to 1/T 



2 
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The heat capacity 



The graph of the heat capacity shows a peak. Beyond the 
peak, the graph falls off as 1/T 2 . Both features are 
characteristic of the magnetic energy levels. 



C 




T 



We shall learn later that heat capacities usually increases. So 
having a peak is unusual. For this reason, the peak is called the 
Schottky anomaly. 

Note that when we say "high temperature" here, it is relative 
to the magnetic energy levels. The Schottky anomaly typically 
occurs around 1 K. 
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The entropy. 



Recall that there is also a formula to derive the entropy from 
the partition function. 

S = Nk B In Z + Nk B T dln Z 



dT 



Substituting the partition function 



Z = 2 cosh(n B B/k B T) 
we get the entropy for the spin-1/2 salt: 



S = Nk B \n 



2 cosh 



-^tanhf^ 
T \k B T 
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The entropy. 



The graph for the entropy looks like this: 



Nk B ln 2 




Again, we would be interested in the high and low temperature 
behaviours. 



Recall the formula for entropy: S = k B In Q 
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Behaviour of the entropy 



At low temperature, most atoms would be in the ground state. 

+HbB 

- K i B B —tf— 

There is only 1 way to arrange this. So the entropy would tend 
to S = k B \n 1 = 0. 



At high temperature, the difference between the energy levels 
would be small compared with the k B T in exp(-e/k B T). 

+H B B — f 

-(ll b B — \ 

The atom is equally likely to be in either level. There are 2 
possible arrangements for each of the TV atoms - 2^ in total. 



So the entropy would tend to S = k B In 2 iV = Nk B In 2. 

PHYS393/Hock 35 2009-2010 



3.2 Simple Harmonic Oscillator 
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The Simple Harmonic Oscillator 



The Simple Harmonic Oscillator is one of the early attempts to 
understand the heat capacity of solids using quantum 
mechanics. 

It was proposed by Einstein. The idea is to think of - to model 
- a solid as atoms vibrating independently of one another other. 



Glazer and Wark (2001) 




Each atom sits in its own potential well. It could undergo 
simple harmonic oscillation in all 3 directions x, y and z. 
Vibration in any one direction is independent of vibration in the 
other two directions. 
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The 1-D case 



From quantum mechanics, we know that the energy for a 1-D 
simple harmonic oscillator is 

£n = (n+ l/2)hu. 

where n is an integer, and v is the frequency. 



Potential energy 
of form 




E n = (n + J-}fiw 



I rite mud ear separation * 



http : / /hyperphysics . phy-astr . gsu . edu/hbase/ quantum/hose . html 
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The 1-D case 



The energy 

s n = (n+ 1/2) hi/, 
can be substituted into the partition function 

Z = $>xp(-e < /fc B T) 

i 

This gives 

zid = E ex p(-( n + 1 / 2 )W^) 

n 

_ exp(-hv/2k B T) 
~ 1 - exp(-hv/k B T) 

The energy can then be found using this formula: 

U = Nk B T* dlnZ ™ 

B dT 
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The 1-D case 



Using the partition function, we obtain the formula for the 
energy: 



2 exp(hv/k B T) - 1 

We are often interested in materials at high and low 
temperatures, because the behaviour can get simpler. 

To obtain the high temperature limit for the energy, note that 
the exponential function: 

e x -+ 1 + x 

for small x. For for high T, 



Nhu 



+ 



Nhu 



U = 




Substituting into the above formula for U gives 



U = 



Nhv 



+ Nk B T 



2 
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The 1-D case 



We have obtained 



Nhu 

U = h Nk B T 

2 



The first term on the right is a constant. The second term 
increases with temperature T. So for high T, the first term can 
be neglected, and we get the high temperature limit for the 
energy: 

U Nk B T as T ^ oo. 



PHYS393/Hock 



41 



2009-2010 



The 1-D case 



The quantity that can be directly measured in an experiment is 
the heat capacity. This can be obtained by differentiating U: 

c = dU 



dT 

Differentiating the energy 

Nhv Nhv 

u = —— + 



2 exp(hv/k B T) - 1 
with respect to temperature gives the heat capacity: 

/9\ 2 exp(e/T) 

C-Nk B ^-j (exp(0/T) _ 1)2 
where 6 — hu/k B . 

We now need to extend this to a real material, which is 3-D 
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The 3-D case 



In the 3-D case, the motion in x, y and z are assumed to be 
independent of one another. 

This means that the energy can be obtained by simply adding 
up the energy in each direction, or multiplying the 1-D energy 
by 3. Likewise for the heat capacity. 

The 1-D heat capacity is: 



where 0 = hu/k B . 

So we simply insert a factor of 3 and get the 3-D heat capacity: 





PHYS393/Hock 



43 



2009-2010 



The heat capacity 



The 3-D heat capacity is: 



C3 D = 3« B |^ 2 * XPWT) 



TJ (exp(0/T) - l) 2 

We would be interested in the high temperature behaviour. 
Using the previous method, when T is large, 

exp J « 1 + - 
The above formula becomes 

C 3D = 3A^) 2 ^ = 3A^(l + i 



For large T, 0/T in the bracket would be a lot smaller than 1, so 

C 3D SA^/c^ as T 00. 
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The heat capacity 



Using the formula we have obtained: 

2 exp(0/T) 



C 3D = 3Nk B (- 



7 (exp(0/T) - l) 2 

we can skectch the graph for the heat capacity. It looks like 
this: 



3Nk 




As we have just shown, the heat capacity approaches 3Nk B at 
high T. 



This agrees with the empirical Law of Dulong and Petit, who 
measured many solids at room temperature in 1819. 
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Compare with experiment 



It also agrees with experiment at lower temperature, 
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A. Einstein, Ann. Physik, vol. 22, p. 186 (1907) 



However, when the temperature gets too low, it falls off 
exponentially. 

This does not agree with experiments, which show that C oc T 3 . 

This is because the effects of co-ordinated movements among 
atoms become important. A new model using phonons is 
needed. We shall learn about this later. 
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The equipartition theorem 



We have seen that at high temperature, the heat capacity 
approaches 3Nk B . 

This can be understood this using the "equipartition theorem," 
which states that: 



"The mean value of each independent quadratic term in the 
energy is equal to \k B T." 

In our case, "quadratic term" refers to the kinetic energy 
l>mv 2 , or the potential energy of the harmonic oscillator, ^kx 2 . 

We can prove this by finding, for example, the mean kinetic 
energy in the Boltzmann distribution: 



f -kmv 2 exp( — irm; 2 lkuT)dv 1 
L2 ^_2_ 1 * = -k B T 

J exp( — j>mv 2 /k B T)dv 
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Explaining the Dulong and Petit Law 



The equipartition theorem can help us the understand the 
Dulong and Petit Law - that the heat capacity is 3Nk B at high 
temperature. 

There are TV atoms. 

Each atom vibrates in 3 directions. 

In each direction, there is kinetic energy and potential energy. 
Each type of energy has a mean value of \k B T. 
The total is [/ = iV x 3 x 2 x \k B T = 3Nk B T. 
Therefore, the heat capacity is C = U/T = 3Nk B . 
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What we have learnt so far 



1. The distinguishable particles model of the paramagnetic ions 
in a solid gives a good description of the magnetisation. 

2. The partition function is given by Z = exp(— ei/k B T). 
This can be used to derive quantities like energy, heat capacity, 
entropy, ... 

3. We can often explain the behaviour of these quantities at 
high or low temperature physically, in a simple way. 

4. The simple harmonic oscillator model of a solid gives a good 
description of the heat capacity, except at very low 
temperature. 
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Some exercises 



3.3 Exercises 
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Some exercises 



Exercise 1 

Given that the total energy of a system is 

U = AJ2^ exp(-s,//e B T) 

i 

where A is defined by the normalisation from 

N = Aj2exp(-e l /k B T) 

i 

demonstrate that 

9 d\nZ 

U = NkuT 

* dT 

where 



Z = 5>xp(-s,//c B T) 
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Some exercises 



Answer 



First, note that since Z = Z)exp(— e^/fe^T) 



TV = A$>xp(- e< /fc B T) = AZ 



Next, 



<9lnZ ldZ 1 ^ Si , n ^ 

= = - Y — l -^exp(-Si/k B T) 

8T ZdT Z^kuT 2 K 11 B J 



B 

Multiplying the two equations above, 



dlnZ A „ , /7 _ N 
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Some exercises 



Rearranging, 



Nk B T 2 ^^ = A^exp(-s^T) 



which is U. Therefore 

B dT 
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Some exercises 



Exercise 2 

The total energy of a paramagnet is 

U = -Nii B Btznh(ii B B/k B T) 

Show that the heat capacity is given by 

C-Nk B {-) (exp((?/r) + 1)2 
where 0 = 2n B B/k B 

Show that the limiting form at high temperature is 
and at low temperature is 

e\ 2 



C = Nk B [^ exp(-0/T) 
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Some exercises 



Answer 2 



The heat capacity is obtained by differentiating 

U = -Nn B Btan\\(iJ, B B/k B T): 



cosh ( 

\k B T 



i -2 



In terms of 0 = 2^ B B/k Bl 



C = 



Nk B 
Nk B 



cosh ( — 



= Nk B 



= Nk 



B 



-) 

2TJ 
(zr) 

») 2 

TJ (exp(0/2T) + exp(-6>/2T)) 2 'exp((9/T) 
0\ 2 exp(0/T) 
TJ (exp(6>/T) + l) 2 



(exp(0/2T) + exp(-6>/2T)) 2 

1 exp(0/T) 
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Some exercises 



At high temperature, exp(6>/T) — »■ 1 

6\ 2 exp(6>/T) 



C - Nk B ^-j (exQ(e/T) + 1)2 



/ 6> \ 2 1 



Vt7 (l + i) 2 

At low temperature, exp(#/T) + 1 — ► exp(#/T) 



\TJ (exp(0/T) + l) 2 



6>\ 2 exp(0/T) 



.77 exp(6>/T) 2 
= Nk B (^J exp(-0/T) 



PHYS393/Hock 



56 



2009-2010 



